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Abstract 

We study a recently proposed gauge invariant, non-local pure gauge observable, 
which is equivalent to the gluon propagator in a certain gauge. The correlator 
describes a gluon coupled to static sources and decays with eigenvalues of the 
Hamiltonian, permitting a non-perturbative definition of a gluonic parton mass. 
Detailed numerical tests of the observable are performed in SU(2) gauge theories 
in 2+1 dimensions. In a Higgs model it reproduces the physical W-boson mass, 
while in a confining theory its non-local nature results in almost exclusive pro¬ 
jection onto torelonic states. However, the gluon mass can also be related to the 
mass difference between the lowest gluelumps, i.e. gluon configurations bound 
to adjoint sources. Its value is computed in SU(2) pure gauge theory in 2+1 
dimensions to be mA = 0.360(19)g 2 , which plays an important role as “magnetic 
mass” of the four-dimensional theory at high temperatures. The same quantity 
is found to determine the mass splitting between vector and scalar mesons in 
the three-dimensional theory. In SU(3) pure gauge theory in 3+1 dimensions 
the corresponding gluelump splitting yields niA ~ 370 MeV. Possible relations 
of this quantity with the QCD heavy meson spectrum are discussed. 
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1. Introduction 


The physics of confinement of quarks in QCD is well established experimentally and 
by numerical simulations of lattice gauge theories. However, little is known about 
the confiment mechanism, i.e. how parton physics at short distances smoothly evolves 
into hadron physics at large distances. Understanding this transition requires non- 
perturbative knowledge of the parton dynamics at all length scales. Such knowledge 
is indispensable in the context of finite temperature and baryon density, where matter 
is believed to appear in a deconfined state whose collective physical properties are 
determined by parton dynamics. The latter is encoded in the Green functions of quark 
and gluon fields, which in general are not gauge invariant. This raises the conceptual 
problem of how to define non-perturbative partons and extract physical information 
about them. In this work an attempt is made to answer this question for a suitably 
defined gluon correlation function. Once this most elementary case is understood, it 
may be generalized to quarks and more complex Green functions. 

In perturbation theory one fixes a gauge and studies partons directly. Although 
field propagators are gauge dependent and not physical observables, they carry gauge 
invariant information about the parton dynamics in their singularity structure. Pole 
masses defined from the gluon and quark propagators have been proved to be gauge 
independent to every finite order in perturbation theory 0.0. However, perturbation 
theory is limited by the requirement of weak coupling. On the other hand, numerical 
results obtained by fixing a gauge on the lattice @ have often been inconclusive or 
controversial in the past. This is due to the difficulty to fix a gauge uniquely and 
avoid the problem of Gribov copies 0. Moreover, most complete gauge fixings (e.g. 
the Landau gauge) violate the positivity of the transfer matrix, thus obstructing a 
quantum mechanical interpretation of the results. Finally, gauge independence of any 
result is not evident, but has to be demonstrated numerically by comparing different 
gauges, while it is extremely difficult to control the mentioned problems for each gauge. 
An overview with references to numerical work may be found in [|5[. 

A non-perturbative argument for a non-zero dynamical mass scale associated with 
gluons was given a long time ago based on the static potential of adjoint sources. 
If gluons were non-perturbatively massless, they could be pair produced at no energy 
cost and the adjoint potential would be screened even for small separations. In lattice 
simulations, by contrast, one observes a linear rise up to rather large distances, until 
the adjoint string breaks at ~ 10 Mq 1 , in units of the lightest glueball |7|, ||. 

Recently, non-local, gauge invariant observables have been constructed in lattice 
gauge theory which are equivalent to the gluon propagator in certain gauges. The 
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transfer matrix formalism has been used to prove that these gluon correlators decay 
exponentially with eigenvalues of the Hamiltonian [|| , implying a pole structure in mo¬ 
mentum space. The energy gap between the vacuum and the lowest excited eigenvalue 
then represents a gauge invariant definition for a gluonic parton mass. 

The purpose of the current work is to compute the value of the gluon mass in pure 
gauge theory and to establish its relation to physical observables. The main results 
are as follows: A gauge invariant probe of the gluon energy is obtained by coupling 
it to external static sources and measuring the energy of the composite object. The 
problem is then reduced to separate the contribution of the gluon field from that of 
the sources. In principle, this can be achieved in two ways. First, by dressing the field 
variables with a non-local functional representing an average over all positions of the 
sources. Second, by introducing explicit fields for the sources, which can be integrated 
over analytically. In both cases the resulting gauge invariant observables can be related 
to the gluon propagator in a particular gauge. The second procedure has the advantage 
of being local and numerically feasible in a confining theory. Moreover, it lends itself to 
a straightforward interpretation of the gluon mass, whose value equals that of the level 
splitting between the lightest “gluelumps”, i.e. gluonic configurations coupled to an 
adjoint source. The latter can be computed straightforwardly without need for gauge 
fixing or introduction of non-local fields. We calculate its value for SU(2) Yang-Mills 
theory in 2+1 dimensions and find rriA = 0.36(2)g 2 . This result plays an important 
role for the four-dimensional high temperature physics, where it corresponds to the 
magnetic mass with g 2 ~ g\ d T. For SU(3) in 3+1 dimensions, the same splitting has 
been computed in the literature to be ttia = 368(7) MeV ]10). In three dimensions one 


furthermore finds the splitting of heavy vector and scalar mesons to be given by the 
gluon mass. The situation in QCD might be similar, but a careful study of the static 
limit is necessary to settle this question. 

The outline of the paper is as follows. In Sec. |2| we review the Hamiltonian formalism 
of lattice gauge theory, which in its strong coupling limit permits a natural definition 
of a gauge invariant mass gap in terms of a gluon coupled to external sources. Sec. [8] 
recalls the main results of p|, relating the Hilbert space definition to a non-local Eu¬ 
clidean observable amenable to numerical simulations. This operator is tested in 2+1 
dimensional theories in Sec. £|. The relation to the gluelump spectrum is derived in 
Sec. ||, and numerical results are quoted for SU(2) in 2+1 and SU(3) in 3+1 dimen¬ 
sions. Finally, Sec. |6] discusses the connection to the heavy meson spectrum for SU(2) 
in 2+1 dimensions as well as QCD, before conclusions are given in Sec. |7. A continuum 
formulation of gluons probed by sources as well as lattice actions and parameters used 
in the simulations are given in two appendices. 
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2. The gluon in the Hamiltonian formalism 


To begin, it is instructive to recall the Hamiltonian formulation of lattice gauge theory 
111 , where it is easy to see that gauge invariant information is associated with held 


variables and how it can be interpreted physically. We are interested in SU(N) gauge 
theory on an L d x N t lattice with periodic boundary conditions. In a Hilbert space 
formulation |TT], [L2], [13| a spatial sublattice L d at a fixed time is considered, with link 
variables f/,(x) and the corresponding held operators C/j(x). The wave functions form a 
Hilbert space Ho of all complex, square integrable functions if)[U] defined on the gauge 
group G: Ho = [L 2 (G)] dLd . A unitary operator R(g) imposes gauge transformations 
g(x) E G on wave funtions according to 


(&/>) [U] = ^[U 9 ], U 9 (x) = g(x)U i (x)g t (x + i). (2.1) 

Wave functions of physical particle states are gauge invariant, iIj[U 9 ] = ?/; [U ], thus 
forming a subspace H cHq- Any wave function 0[ U ] E Ho can be projected on the 
physical subspace by means of the projection operator P, 


(Ptl>){U] = Tldg(x) i>{U°]. (2.2) 

J X 

The Kogut-Susskind Hamiltonian acting on the full space Hq is obtained by quan¬ 


tizing the theory in temporal gauge U 0 — 1 and reads [XT 

g 2 2 


^ £W(x )) 2 - E ReTr P„(x), 
^ - y a *,i<j 


(2,3) 


where f),(x) — t/,(x)r/ J (x + i)[/'/(xH-i+j)f/j(x) is the elementary plaquette held. The 
components of electric held and link operators satisfy the commutation relations 


£?(x),£?(y) 

A a (x),tW y) 


i/abc-%( X )<M xy , 

-T a U i , af} (x)5 ij 8 xy . 


(2.4) 


In principle the whole configuration space may be constructed by applying link opera¬ 
tors in all possible combinations to the vacuum state, 0 0 = 1. In particular, doing this 
just once we have 

^apl u } = Ui, a g(x)^ 0 = £Wx)0 o . (2.5) 

Thus, 0 1 [17] is the wave function corresponding to one link variable being excited 
and all others being in the vacuum configuration. Clearly, it transforms non-trivially, 
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^a^ 9 ] = da'y (x)< 7 £ 5 (x + i )if)js[U\- Gauge invariant wave functions are obtained by 
exciting closed loops of links, the simplest being a plaquette, A P [U} = P^-(x). 

In the limit of strong coupling, a^ 3 ~ d ^g 2 —> oo, the potential term in the Hamiltonian 
decouples. Although this limit is far from the continuum and the physical situation, it 
illustrates the structure of the Hilbert space Rq. In the strong coupling limit A 1 is an 
eigenstate of the Hamiltonian, 

itoKeVA = E 1 = E^C f , C f — (TV 2 — 1)/2N. (2.6) 

Excitation of other links increases the energy by one unit for every link. The wave 
function for a string of electric flux between charges at x, y = x + Zi, 

€p[U] = [G 8 (x)G,(x + i)... I/ 8 (x + (/ - l)i)U = Mx, y), (2.7) 

is an eigenstate with energy IE \, whereas 'ip P [U] has energy 4ifj. The plaquette wave 
function is in the Hilbert space of gauge invariant functions, ip p [U\ G 7i. On the other 
hand, the wave functions A 1,5 are non-invariant and hence not elements of the projected 
Hilbert space. Nevertheless, because time independent gauge transformations commute 
with the Hamiltonian, [R, H] — 0, the eigenvalues belonging to these wave functions 
are gauge invariant. 

It is now natural to define a “gluon mass gap” by the energy it takes to excite one 
gluon configuration out of the vacuum. Defining the usual lattice gluon held 


i 

I/Ax) - U\ (x) - -J-Tr (l 

J Ax) — U- (x)) 

2 g 

n ’ 1 v ’ N V 

f' \ / l \ / J 


one has for the corresponding wave function ip A s[U] = A,a/3( x )Ao hi the strong coupling 
limit 

tfo tfelU] = E^{U]. (2,9) 

Hence, in the Hamiltonian strong coupling limit the gluon has a gauge invariant mass 
gap, which is an eigenvalue of H 0 to a gauge covariant eigenfunction. Note that the 
wave functions for a unit of flux and the gluon wave function are constructed from 
the same link variable. The difference is in the transformation property. i!j a picks 
the traceless part of the link variable only, and hence is an (iV 2 — l)-plet of the form 
~ c a [U]T a with some real functions c a [U]. 

Away from the strong coupling limit the energy levels and the precise form of the 
eigenstates will change. However, their transformation properties remain the same since 
gauge transformations commute with the full Hamiltonian. Quite generally, eigenstates 
with non-trivial transformation behaviour describe strings of flux with energy E. They 
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can be made manifestly gauge invariant and “physical” by letting the flux end on static 
sources. The total wave functions then are gauge invariant combinations of the wave 
function for the source, 0 ff (x) = g(x)0(x), and the gauge part, 

*f[e] = C(x)^(x + w'Jlu], Klu] = Y g (API - *iPl) ■ (2.10) 

In summary, in addition to the gauge invariant particle states of 7d, 7do contains 
additional states describing static potentials, or gauge held configurations probed by 
external charges. The corresponding pure gauge wave functions transform non-trivially, 
but the associated energies are fully gauge invariant and contain valuable information 
about the dynamics of charged states. In this respect a gluon wave function is on the 
same footing as a wave function for a flux element, whose traceless part it constitutes. 
We conclude that a gluon mass may be defined by the minimal energy it takes to 
excite a flux element in an (iV 2 — 1) or adjoint state. This energy is not related to an 
asymptotic particle state, but rather a held coupled to static sources. The expression 
Eq. (|2.1Q ) suggests that it is a short distance property of the static potential with the 
sources and the hux in an adjoint state. 


2.1. Transfer matrix and Euclidean correlation functions 


In order to make the spectrum amenable to numerical simulations, the Hilbert space 
picture has to be connected to Euclidean correlation functions by means of the transfer 


matrix | p~2| . |T3f 

time, 


The latter is an integral operator translating the wave functions in 


i>[u <+i] 


( tw ) [Ut] 

J n«..(x) exp— St[Ut+uUt]^[U t ], 


( 2 . 11 ) 


where S t is the Wilson action of two neighbouring timeslices in temporal gauge. T 0 is a 
bounded, self-adjoint operator with a strictly positive spectrum JE|, allowing to define 
a bounded Hamiltonian Hq = — 1/alnTo, which up to corrections 0(a 2 ), is identical 
to the Kogut-Susskind Hamiltonian Eq. Q2.3|). Through the projection 


T = PTn 


( 2 . 12 ) 


one defines a transfer matrix T acting on the gauge invariant subspace only, as well as 
the corresponding Hamiltonian H. 

Writing the quantum mechanical trace over a complete set of states on the space 
Tt as Tr O = J2 n ( n \0\ n )i a one t° one relation between Euclidean expectation values 
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and quantum mechanical traces can be established. For correlators of gauge invariant, 
local operators 0(x) the result is 


(C?(x,O)0(x,t)) = Z~ lr fr {T N ~ t O(x.)T t O(x)} 

Nj —^° l(0|C , (x)|n)| 2 exp ~(E n - E 0 )t . (2.13) 

n 


The eigenvalues and states |n > are those of the Hamiltonian H and describe the 
asymptotic particle states, in the pure gauge theory the glueballs. 

On the other hand, the static potential is probed by the Wilson loop, with the 
temporal lines representing the external charges. In this case the correspondence is 

<W(|x — y|, f)) = Z- 1 tr{f»-‘f)„ () (x,y)T' 0 ‘(7tyx,y)} 

y |{0|(7 oS (x,y)|n)| 2 exp-(B„ - E„)t . (2.14) 


Here the |n > and E n are eigenstates and eigenvalues of Ho and do not correspond to 
particle states, but rather to field energies in the presence of external charges. Quite 
generally, correlators involving temporal Wilson lines describe configurations involving 
static charges, their expontential fall-off is dictated by Ho, and the intermediate states 
transform non-trivially. 

Euclidean and Hilbert space formalism are but a different language for the same 
physics. Following the observations in the last section, it must then be possible to con¬ 
struct a Euclidean operator 0[U] whose expectation value can be expressed analogous 
to Eq. fl2.14p with f/(x, y) —» A(x). 


3. A non-local, gauge invariant gluon operator 

It has recently been shown that such a Euclidean expression can be realized with an 
auxiliary complex TV-plet f[U ] transforming in the fundamental representation of the 
group || . This procedure is not unique. One possibility is to use eigenfunctions of the 
covariant Laplacian, which is a hermitian operator with a strictly positive spectrum, 

-A = £ [2/l’lW - (i + A) - 0j(. - - A) 

n 

= A n f {n \x), X n > 0. (3.1) 

Its eigenvectors have the desired transformation property f^ 9 (x) = g(x)f( n \x). They 
provide a unique mapping U —> f[U} except when eigenvalues are degenerate or |/| = 0. 
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In practical simulations the probability of generating such configurations is essentially 
zero [[TJ- These properties have been used previously for gauge fixing without Gribov 
copies ||TI| and to construct blockspins for the derivation of effective theories [[HJ . The 
eigenvectors are non-local in the sense that they depend on all link variables. In order 
to maintain the transfer matrix formalism the f(x) have to be local in time. This is 
achieved by considering the spatial Laplacian A 2 \Uj\ in Eq. (|3.1|) , which then is defined 
in every timeslice individually and independent of Uq. 

The eigenvectors are used to construct an N x N matrix fi(a;) G SU(N) following 
Since Eq. (|3 . 1|) only determines them up to a phase, this leaves a remaining 
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freedom in Q(x). In the case of SU( 2), to every eigenvector / there is a degenerate 
second one given by its charge conjugate ir 2 /*. In order to have a smooth held the 
solution to the lowest eigenvalue is chosen in practice. These may now be combined 
into 


£)(#) = 


/i (0) W 


l/ (0) MI l /f M 


f (0)* 
.12 

fl 0> ( 


X 


X 


G SU (2). 


(3.2) 


The whole matrix satisfies the Laplace equation, —Af[U]£l = A°fi, and so does f 1h, 
where h may be any global SU( 2) matrix. For SU(3) there is no degeneracy of the 
eigenvalues in general. In this case one solves for the three lowest eigenvectors to 
construct the matrix 12 = which is then determined up to a factor 

h = diag(exp(io;i), exp(ico> 2 ), exp(ic<j 3 )), Xj* = 0. (For an alternative construction and 
a numerical implementation, see [I7j]). The remaining indeterminacy h may be different 
for every timeslice. This is summarized by the transformation behaviour 


n 9 (x) = g(x)fl(x)h\t), 


(3.3) 


where h(t) is free. 

We can now define composite link and gluon fields 


V,(x) 

A^(x) 


Vt\x)U^(x)Vt(x + fi), 


i 

2g l 


V F (x) - - -IV (U(i) - V}(x)) 


(3.4) 


both transforming as Of(x) = h(t)Oi{x)h)(t), whereas Vq(x ) = h(t)Vo(x)hJ(t + 1). 
Hence the Ai(x) are gauge invariant under spatial transformations g(x), but transform 
under time-dependent rotations corresponding to the residual symmetry of the spatial 
Laplacian. Note that the functions H and hence the composite gauge fields have no 
definite Zjv symmetry and couple to all sectors. 

With the abbreviation 

< 2-1 

V 0 (x;t 1 ,t 2 ) = n HOMO (3.5) 

t'=ti 
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for a temporal Wilson line, we can now construct the gauge invariant Euclidean oper¬ 
ator 


0[U] = Tr [A(x,0)V'o(y;0, t)A t (x,t)V^(z;0,t) 


(3.6) 


It represents a correlator in t of the composite field Ai, where Wilson lines are inserted 
to ensure full gauge invariance. Note however, that these may be placed at any y, z. For 
the particular choice x = y = z, the Wilson lines merge into one adjoint representation 
line V^ ab (x; 0, t) = Tr [T a Vb(x; 0, t)T b V^ (x; 0, t)], 


0[U] = -A a (x,0)E 0 A afe (x;0,t)4 b (x,t). 


(3.7) 


Eq. (|3T7|) represents a correlator describing a gluon bound to an adjoint source, in 
analogy to the “gluelump” operators JT8[. In Appendix [X] the corresponding continuum 
expression is derived by coupling the gluon field to explicit fields for the static sources, 
which are then integrated out analytically. 

In H the transfer matrix formalism has been employed to convert Eq. (|3.6|) into a 
quantum mechanical trace, 


(< 0[U]) = Z -'Tr {ff‘„ ii,^(x)} , 


(3,8) 


where we have defined a modified transfer matrix Tlo by the ‘Laplacian temporal’ 
gauge in which Vo(x) = 1. It has been proved that Tlo has the same spectrum as T 0 
[:9|. Initially we had only known that T 0 is invariant under time independent gauge 
transformations. Now we have the stronger result that T 0 and Tlo have the same 
spectrum, even though they are related by a time dependent gauge transformation. 
This finding is the lattice analogue to a continuum quantization of the Schrodinger 
functional: equivalent results are obtained in temporal and any other gauge that is 
invariant under spatial rotations [T^ . 

Hence the expectation value of our operator decays exponentially with eigenvalues 
of the Hamiltonian Hq, 


lim (0[U]) = |(0|A iQ/ g(x)|n)| 2 exp-(E n - E 0 )t 

N t >oc „ 


(3,9) 


Eqs. 


and (|3.9|) are the desired gluon analogue to Eq. Q2.14| ). 


3.1. Existence of a continuum limit 

When the continuum limit is approached, the exponents extracted from Eq. Q3.9|) di¬ 
verge because of the self-energy contributions of the temporal Wilson lines, Fig. [lj. In 















p-q 


q 


p 


Figure 1: The Wilson line self-energy. Eq. ( j 3.1 Of enforces p = p — q = 0. 


order to retain a finite continuum limit, the operator has to be modified such that no 
divergent mass renormalization is present. 

This can be achieved by observing that the transformation behaviour of Vq (x, t) in 
Eq. ( |3.4| ) is independent of the spatial coordinates. In the construction of the operator 
Eq. (|3.6| ), instead of Vq we may then use its timeslice average 


?oM = £v r o(x,<)/ll£v'„(x,oil, 

X X 


(3.10) 


which has been projected back into the group. The timeslice average corresponds to 
the Fourier transform of Vo with zero momentum. If this is done in every timeslice, 
the sources represented by the Vo cannot emit a gluon at one t and reabsorb it at some 
later t as in Fig.|l]. Hence, the mass renormalization of the static source is switched off, 
and the sources remain classical external fields. The presence of fields Vo then merely 
accounts for the transformation behaviour, but has no effect on the gauge held energies 
measured by the modified operator 


(0[U)) 



A(x, 0)Vo(0, t)A(x, t)V$ (0, t) 




(3.11) 


which has a spectral decomposition as in Eq. ( |3.9| ). The energies extracted from the 
expectation values of Eqs. ( |3.6| ) and Q3.11p should then differ by a cut-off dependent 
shift due to the selfenergy contribution Fig. |T|. A non-perturbative parton mass for the 
gluon may now be defined by the gap between the lowest excitation energy and the 
vacuum, 

rrtA = E\ — Eq. (3-12) 

By Fourier transformation of the spectral representation Eq. (|3.9|) it follows trivially 
that the energy eigenvalues appear as poles in momentum space. It depends on the 
behaviour of the higher energy levels whether in the infinite volume and contiuum 
limits the lowest eigenvalue remains an isolated pole or turns into a cut. 
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3.2. Non-uniqueness of fl and the relation to gauge fixing 


How do the above results depend on the particular construction of H? Clearly, any 
0[U] G SU(N) local in time and satisfying the transformation law Eq. (|3.3|) , permits 
construction of a gauge invariant observable Eq. (|3.11| ) . From the spectral representa¬ 
tion it follows that all such observables fall off with the same spectrum, while 12 only 
enters the matrix elements and thus influences the overlap of the operator with the 
eigenstates. 

Of course, the construction of the composite link variable Eq. (!T1) may also be 
viewed as fixing Laplacian gauge on each timeslice [15 . This gauge is incomplete, 
with a global factor h remaining unfixed between time-slices. It can be completed by 
imposing the further condition Vo (f) = 1, thus fixing h(t). In this particular gauge 
the operator Eq. ( |3.11| ) is equivalent to a gauge fixed gluon propagator, falling off 
exponentially with eigenvalues of the transfer matrix. Since the spectrum is unaffected 
by the particular construction of 12, this statement holds for all gauges employing a 
unique Q [U ] local in time. For example, the standard Coulomb gauge is defined by an 
12+) that minimizes the functional 


m = e 


1 -jv Ttw 


fit, 


X ) 


(3.13) 


in every timeslice. 12+) has the desired properties and a residual freedom h(t) G 
SU(N), to be fixed in the same way. A recent implementation is in |^D|. (Of course, 
this gauge condition has the problem that it does not determine 12 uniquely ||). Hence, 
any gauge fixing done in this manner is equivalent to coupling the flux to sources and 
averaging over all their positions. In momentum space this non-perturbative result 
translates into what is also found in perturbation theory: the pole of a propagator is 
gauge independent, whereas the residue, viz. the matrix element, is not. On the other 
hand, Landau gauge is non-local in time, no positive transfer matrix is defined and 
non-perturbatively it is not guaranteed that it probes the same spectrum. 


4. Numerical tests for SU(2) in 2+1 dimensions 

This section is devoted to numerically test the operator constructed in the last section. 
This first exploratory study is done for SU(2) in 2+1 dimensions, for its significantly 
lower numerical cost and fast continuum approach permitting to extract conclusive 
continuum results. In this case the coupling constant g 2 has dimension of mass and 
provides the scale in which all physical results are expressed. 
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4.1. W-boson in a Higgs model 


It is useful to first test the operator Eq. ( |3.11|) in a situation where the result is known 
and controlled by perturbation theory. We therefore consider a SU(2) Higgs model with 
scalar 2x2 matrix fields <h(a;) in the fundamental representation and the continuum 
action 


5 = d s x Tr 


~F l3 F l3 + (A<h) t A$ + m 2 0 &$ + 2A($ t $)' 


Its physical properties are determined by the two dimensionless parameters 

A ml 

x = ~2> y = — • 

r g 4 


(4.1) 


(4.2) 


The lattice action and parameters as well as their relation to the continuum ones are 
given in Appendix |B[. 



Figure 2: Schematic phase diagram of the Higgs model Eq. 0. Pure gauge theory is 
reached in the limit x,y —> oo. 


The phase diagram of the theory has been determined non-perturbatively in the 
continuum limit and is shown schematically in Fig. [|: Higgs and confinement regions 
are separated by a line of first order phase transitions, which ends in a critical point . 
The spectrum in the respective regions as well as its continuous connection through 
the crossover region has been studied extensively [^, [23| . 


In the Higgs region, the theory has a physical vector boson carrying the quantum 
numbers of the gluon, coupling to the gauge invariant composite operator 


Wf(x) = Tr ^T a< h' l ’(a;)iAj < f)(a;)) (continuum), 

W?(x) = Tr (T a &(x)Ui(x)<S>(x + i)) (lattice). (4.3) 
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Its mass is determined by the decay of the corresponding correlation function 

(W i “(x,0)W i Q (x,f))~e- M ^. (4.4) 

My is an eigenvalue of the Hamiltonian acting on the gauge invariant subspace 7i, 
and the W-boson is an asymptotic particle state. On the other hand, in perturbation 
theory one fixes a gauge and computes the vector boson mass from the fall-off of the 
gauge field propagator, 

(4,5) 

In momentum space the mass corresponds to the renormalized pole of the propagator. 
At a finite order in perturbation theory one finds rriA = My. With our operator 
constructed in the last section, Eq. ( |3.11|) , we now have a lattice implementable gauge 
field correlator without any reference to the scalar fields. We non-perturbatively know 
it to decay exponentially, where is an eigenvalue of the Hamiltonian acting on the 
entire Hilbert space Ho- 

Fig. |3] shows a comparison of the ground state mass as obtained from the two corre¬ 
lation functions, measured at one point in the Higgs phase, x = 0.0239, y = —0.02, for 
different lattice spacings. Lattice sizes were chosen large enough for the masses to be 
free of finite volume effects, based on the results of [^]. Within statistical errors, the 
two operators yield identical results. The mass extracted from the new operator indeed 
extrapolates to a physical continuum limit, as promised in Sec. [|. As expected for an 
object with dimension mass 2 , the lowest eigenvalue of the spatial covariant Laplacian 
diverges strongly in the continuum limit. For the correlation function this does not 
matter because the eigenvalue itself does not enter the definition of the correlation 
function. 

We then conclude that the gauge field propagator has a non-perturbative pole in 
the Higgs region of the model, and within our accuracy we find mn = My also non- 
perturbatively. While this may not be surprising, we stress that it is a non-trivial 
result. What has been known is that the perturbatively defined pole in the field 
propagator gives a good approximation of the physical W-mass in suitable gauges 
and when the perturbative Higgs expectation value is large compared to quantum 
fluctuations [|24], |25j, i.e. when perturbation theory is reliable. However, it has not 
been clear whether this pole also exists non-perturbatively, since that requires a non- 
perturbative definition for Moreover, it is important to realize that, because of 
the residual gauge freedom in the composite fields, Eq. (|3.4|) , we have 

(A‘[U]W‘{4>,U]) = 0, (4.6) 

i.e. the operators are orthogonal, and so are the corresponding eigenstates of the Hamil¬ 
tonian. This is yet another reflection of the fact that the eigenstates of the Hamiltonian 
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Figure 3: Left: The W-boson mass in a 3d Higgs phase, computed from the gauge invari¬ 
ant composite operator Wf[4>, U], Eq. (\j.3f ), and the gauge field propagator, Eq. (\3.11\) . 
The data for My are from jll/ . Right: The lowest eigenvalue of the lattice Laplacian, 
Eq. ([OJ. 


Ho contributing to the correlator Eq. (|3.11|) are in a different sector of the ffilbert space 
than the eigenstates contributing to the correlator Eq. (ED- Consequently, tha, My 
are really different eigenvalues of H 0 which are degenerate in a ffiggs dynamics only. 

The analytic connectedness of the phase diagram implies that there is a one to one 
mapping of the the entire spectrum of the Hamiltonian Ho between Higgs and con¬ 
finement regions of the phase diagram For example, the W-boson state described 
by the operator Eq. ( |4.3| ), being an eigenstate of the projected transfer matrix T, be- 

As we shall see, 


comes a vector meson bound state in the confinement region [^, [ 
rriA 7 ^ My in a confining regime, as is also found in resummed perturbation theory 
|27| and in lattice Landau gauge [|8j. The parameters can be continuously tuned 
to x, y —> oo, where the scalar fields are infinitely heavy and decouple to leave us with a 
pure gauge theory. All meson states including My disappear from the (finite) spectrum 
in this limit. The question then is what happens to the pure gauge quantity rriA- 


4.2. Yang Mills theory 

In addition to a gluon, one may also construct other operators from the composite links 
Eq. (|3.4| ). In particular, simply taking the trace, 

O v (x) = Tr(V-(aj)), (4.7) 
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produces an operator coupling to 0 ++ states. It is fully invariant under gauge trans¬ 
formations, the residual freedom h in the definition of the composite link discussed in 
Sec. |3] drops out under the trace. Hence, the correlation (Oy(x)Oy(y)) falls off with 
the spectrum of the projected Hamiltonian //, and one expects to be able to extract 
the 0 ++ glueball mass from it. 

Fig. [| shows the masses in lattice units extracted from the correlators in the 0 ++ 
and 1 channels, respectively, plotted against the spatial length of the lattice, L. At 
intermediate and large distances we see a linear rise with L, which is reminiscent of 
torelons, i.e. flux loop states winding around the lattice |29|. The well-known 0 ++ 
torclon couples to Polyakov loops P- L \x) = II n=o Ui(x + ni), winding around a spatial 
direction of the lattice. Its mass extracted from the corresponding correlation function 
is related to the string tension in a finite volume, 

^Tr (x, 0) Tr P/ L) (x, t)^ ~ e ~ am P (L)t^ am P (L) = a 2 a L L. (4.8) 


Torelons are physical states in the gauge invariant sector in a finite volume, but be¬ 
come infinitely heavy and decouple when the infinite volume limit is taken. As Fig. ^ 
demonstrates, the two operators Oy, Tr P dd give identical results and we conclude that 
Oy projects onto the torclon. 

This suggests to interpret the lower 1 state as torelonic state as well. In order to 
confirm this, we compare with the traceless part of the Polyakov loop. The latter is not 
gauge invariant, but if we replace the spatial link by a spatial Polykov loop, U t —> P^ 
in Eq. o, the correlator Eq. (|3.11|) is again gauge invariant. The resulting masses are 
also shown in the plot and found to be identical to the ones from the original correlator. 
Hence we conclude that we have found a new eigenstate of the Hamiltonian for a pure 
gauge theory in a finite volume, corresponding to a torelon in an adjoint state. Because 
it transforms non-trivially, it is not an asymptotic state like the standard torelon, but 
rather an eigenstate of Ho. 

For L = 24, we have also calculated the correlators with an il[U] computed from 
the Coulomb gauge fixing condition Eq. ( fh 13|) . As Fig. [| shows, the results are, in 
accordance with our theoretical considerations, equivalent to those obtained from the 
Laplacian procedure. In particular, since the Q emerging in Coulomb gauge is also 
non-local, it projects on the same winding states. The small deviation observed in the 
1 mass might be related to the influence of Gribov copies. 

Fig. [3| is essentially unchanged when computed in the confining regime of the Higgs 
model rather than in the pure gauge theory. This is in agreement with previous findings 
that the confining gauge sector of the Higgs model is almost insensitive to the presence 
of the scalar fields [^, ^3|. It confirms that the gluon correlator Eq. (|3.11|) does not 


project on the vector meson with mass My, but rather on a pure gauge quantity. 
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Figure 4: Effective masses of 1 and 0 ++ operators using the composite link V[U\ and 
Polyakov loops, as described in the text. The simulation is in YM-theory at (3 — 9 with 
spatial lattice size L. 

4.3. Localized vs. torelon states 

Does the adjoint torelon state have anything to do with the gluon? For finite L this 
is not the case, just like the 0 ++ torelon is not related to the glueball. The analytic 
connectedness of the phase diagram holds for the theory in infinite volume. This implies 
that the infinite volume pole in the gluon correlator found in the Higgs regime cannot 
disappear from the spectrum when the parameters are changed to the confining regime. 

The reason for seeing only torelonic states is a projection problem of the maximally 
non-local operators: the functional 0[L/] depends on all links in the spatial volume. In 
a Higgs regime this is of no import because colour is screened and the links far from 
each other do not communicate. In a confining regime on the other hand, adjacent 
links form flux tubes stretching over the whole spatial lattice. An operator constructed 
from D must have a large overlap with such non-local objects, and a very small overlap 
with localized states. Put in another way, in the composite gauge field the link to be 
correlated, U,fx), is “drowned” by the dynamics of all other links contributing through 
fl[£7], a situation that worsens as the spatial volume is increased. For example, at 
some L the 0 ++ torelon becomes heavier than the lightest glueball, which then should 
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dominate the correlation function of Oy. At the lattice spacing chosen in Fig. [|. the 
lightest glueball mass is aM G = 0.755(7) [|3T|. At L = 32 the torelon is heavier, yet 
Fig. f| demonstrates that Oy is blind to the localized glueball state. 

This finding is corroborated by the fact that the problem is absent when correlations 
other than those of pure gauge operators are studied. For example, propagators of an 
adjoint scalar field <&(x) coupled to the gauge fields have been studied in the Coulomb 
gauge in [!2Cf] . No significant finite volume dependence was found for the scalar propaga¬ 


tor. While the corresponding composite operator fT$(:r)Q(:r) also employs a non-local 
function Q[I7], it is local in the scalar field to be correlated. 

We then conclude that pure gauge operators constructed from the composite links 
Vi in a confining dynamics have very weak overlap with localized states, and in their 
current form are not very useful to extract infinite volume physics. Nevertheless, a new 
eigenvalue of the Hamiltonian in the sector with quantum numbers of the gluon has 
been found. Just like the 0 ++ torelon is related to the potential of static charges in 
a singlet state, the adjoint piece of the flux tube should be related to the potential of 
static charges in a triplet state (octet for SU(3)) [02], This question is beyond the 


scope of the present paper and will be pursued elsewhere. However, having extracted 
the energy of a flux tube in an adjoint state, note that it extrapolates to a non-zero 
value ~ 0.32 g 2 for L —> 0. In our Hamiltonian considerations we have interpreted the 
gluon mass as the minimal energy of precisely such a piece of flux. However, one would 
expect this value to be affected by finite volume effects, as one side of the lattice has 
been shrunk to zero. 


The gluon mass and static mesons 


In order to circumvent the projection problem encountered in the previous section, we 
now seek to compute the energy eigenvalue of interest from local observables coupled 
to sources. Rather than averaging over all their positions, the contribution of the 
sources will be computed separately and subtracted. This approach will also provide a 
relation between the gluon propagator in temporal gauge and the ratio of static meson 
correlators, leading to a physical interpretation of the gluon mass in terms of level 
splittings. 

For this purpose we go back to the expressions of a gluon coupled to external sources 
in the continuum, Eq. ( [Ajl ,8| ) from Appendix A, and on the lattice, Eq. (|3.7| ) . For 
a non-perturbative evaluation, we want to avoid the non-local functions used in the 
definition of the composite lattice gauge field. This may be achieved by observing 
that the energies governing the decay of a correlator are completely determined by 
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the Hamiltonian and the quantum numbers of the operators. Instead of discretizing 
the gauge held, we may then equally well employ a local higher dimension operator 
sharing the same quantum numbers, such as the adjoint part of the linear combination 
of plaquettes in the spatial plane transverse to the Wilson line, 


O df (x) = Pij(x) + P-ij(x) + P-x-i(x) + P-j,i{x ) - h.c. 


(5.1) 


In the continuum limit, the operator approaches the covariant derivative of the held 
strength, DjFij(x), and couples to J pc = 1 states. We are then considering the 
two-point function 

(0“(x,0)C/^ t (x;0,()0 4 (x,«)) ~ (5.2) 

with an adjoint temporal Wilson line U ( f ah (x: 0, t). It has the general form of a gluelump 
correlator, describing a bound state of dynamical glue and a static adjoint source ]18| . 
Hence we are looking for the mass of the lightest 1 , or “electric”, gluelump. 

Because we are now working with undressed temporal links rather than with com¬ 
posite ones, the zero momentum projection Eq. ( |3.10|) cannot be applied, and the 
observable still contains the divergent contribution from the adjoint Wilson line and 
its self-energy. We now need to compute this contribution separately and subtract it 
in order to get a finite result in the continuum limit. The Wilson line can be made 
gauge invariant by closing it through the periodic boundary in the t-direction, in which 
case it becomes an adjoint Polyakov loop, whose exponential decay with the temporal 
lattice size N t we need to compute. On the other hand, the Polyakov loop couples to 
adjoint states in all J PC sectors, and asymptotically its exponential decay is the same 
as that of Eq. (|5.2|), where 0(x ) is the operator leading to the smallest coefficient Mo- 


Hence, we are looking for the mass of the lightest gluelump, which is well known |18 
to be the “magnetic” one described by the clover field, 


O f (x) = P id (x) + Pj-i(x) + P-i,-j(x) + P- j}i (x) - h.c. 


(5.3) 


In the continuum, this operator approaches the field strength Fij(x), coupling to 0 
and l +_ states in 2+1 and 3+1 dimensions, respectively. 

The lowest eigenvalue of H 0 determining the asymptotic decay of the gluon correla¬ 
tion function Eq. (|3.11|) should then be the same as the gluelump mass splitting 


TTij\ — M df — M f 


(5.4) 


Note that the gluelump states may be viewed as the zero distance limit of static po¬ 
tentials with the string between the sources being in an excited state specified by 


appropriate quantum numbers [181. This is in accordance with our interpretation of 
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the gluon mass as a short distance property of potentials put forward in the Hamilto¬ 
nian strong coupling analysis, Sec. |2|. The calculation of gluelump splittings involves 
only pure gauge operators local in space and is straightforward to evaluate on a lattice, 
with well controllable infinite volume and continuum limits. 


5.1. Numerical result for SU(2) in 2+1 dimensions 


Gluclumps in three dimensions have been simulated previously in different contexts 
[34 , 35]. Here we are specifically interested in a continuum extrapolation of the mass 


splitting Mdf~M p. The numerical results of a calculation in SU(2) pure gauge theory 
close to the continuum are displayed in Fig. ||. The choice of lattice sizes follows the 
results of [35], where it has been explicitly checked that the lowest gluelump is free of 
finite volume effects. The figure displays the onset of the weak logarithmic divergence 



Figure 5: Masses for the electric (M DF ) and magnetic (M F ) gluelumps, as well as 
their mass splitting. Calculations have been performed for (3 = 12,18, 25 on lattices 
with L = 30,48, 64, respectively. The mass splitting is extrapolated to the continuum. 

in the gluelump masses, as well as the expected scaling for their mass difference. A 
linear extrapolation of the data in 1/(3 to the continuum then gives the final result for 
the SU(2) gluon mass in 2+1 dimensions 

m A = 0.360(19)g 2 , x7 dof = 0.47. (5.5) 
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Comparing with the lightest scalar glueball, Mq = 1.584(17) <7 2 m , one finds Mq /rn A = 
4.4(3). 


5.2. The magnetic mass 


Besides being an interesting test case, the gauge theory in 2+1 dimensions plays an 
important role in four-dimensional physics at finite temperature. Representing the 
Matsubara zero mode sector of the latter, it constitutes the effective theory describing 
all static physics at asymptotically high temperatures. A long standing problem of 
thermal field theories are the severe infrared divergencies encountered in perturbation 
theory. In particular, loops of magnetic gauge fields exhibit divergences on the 
three-dimensional scale g 2 ~ g 2 d T. These may be cured by dynamical generation of a 
“magnetic mass” which, however, is entirely non-perturbative and receives contribu¬ 
tions from all orders in a loop expansion [|36|. Here we have shown that such a magnetic 
mass can indeed be defined non-pertubatively, and that for hot SU(2) gauge theory its 
asymptotic high temperature value is m A = 0.36(2)gf d T. 



ref. 

m A /g 2 

1 -loop gap eq. 

H 

0.38 


27, 38 

0.28 


3S 

0.25 

2 -loop gap eq. 

m 

0.34 

lattice MAG 

42 

0.51(6) 


Table 1: Comparison of magnetic mass values from gap equations and gauge fixed 
lattice simulations. 


In the past, gauge invariant resummation schemes have been designed to compute 
the pole of the gluon propagator self-consistently in three dimensions ||27||, [^7[] -p0| . 


In a Hamiltonian analysis of the three dimensional gauge theory a gauge invariant 
composite gluon variable has been constructed, which in the weak and strong coupling 
limits yields a gluon mass gap as the lowest eigenvalue of the Hamiltonian |4l|]. One 
may now compare the results of these approaches collected in Table |T] with the full 
answer and get an estimate for the quality of the approximations involved. It appears 
that the resummations lead to reasonable answers. In particular the two-loop result 
exhibits convergence towards the result of the last section. Note that lattice simulations 
employing maximal abelian gauge fixing provide numbers that are somewhat off our 
result, underlining the difficulty of obtaining stability in such calculations. I 11 Landau 
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gauge the correlator becomes negative at large distances, and it appears difficult to 


extract a mass 2C . 


5.3. Numerical result for SU(3) in 3+1 dimensions 

The considerations to relate the gluon mass in pure gauge theory to the splitting of 
magnetic and electric gluelumps are directly applicable to SU(3) Yang-Mills theory 
in 3+1 dimensions. The only difference is that the lowest magnetic gluelump and 
the corresponding operator in four dimensions has J PC = l +_ . The mass splitting 
between these lowest states has been computed in [|10|] , which quotes for the continuum 
extrapolation 

niA — M df — M f — 368(7)MeV. (5.6) 


Comparing with the lightest scalar glueball, M G [ 0 ++ ] = 1730(130) MeV |0, we then 
have Mc/rriA = 4.7(5), which is strikingly similar to our 2+1 dimensional calculation in 
SU(2). This is not unexpected, given the close resemblance of other properties of pure 
gauge theories in 2+1 and 3+1 dimensions. Together with the ratio M G /mA , Table 
D collects results for the mass of the lightest scalar glueball and the string breaking 
distance for the flux tube between adjoint sources, suggesting that the confining low 
energy physics is very similar in these theories. 

Similar to the situation in three dimensions, a larger value niA ~ 600 MeV has 
recently been estimated from the gluon propagator in Laplacian gauge [05|. However, 


while this gauge is free from Gribov copies, no positive transfer matrix exists and hence 
the functional form of the decay is not exactly known. 


Table 2: Comparison of some continuum extrapolated pure gauge quantities between 
2+1 and 3+1 dimensions. _ 



2+1 3+1 

SU(2), M G [0 ++ ]/y/a 
SU(3), M g [ 0 ++ }/y/a 
SU(2), M G rt{ ak 
SU(2), M G /m A 
SU(3), M G /m A 

4.72(4) pi 3.74(12) g§ 

4.33(4) 0 3.64(9) 0 

10.3(1.5) 0 ~ 9.7 § 

4.4(3) 

4.7(5) b ni 
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6. Mass splittings between heavy mesons 


In this section we ask what happens in the vicinity of the static limit, when heavy 
but dynamical matter fields are present, and whether the static limit can be attained 
smoothly. One might hope that the gluon mass is then observable in the spectrum of 
sufficiently heavy mesons. Two candidate splittings appear as possibilities. 

First, hybrid meson states are, in the non-relativistic approximation, bound states 
in gluonic excitations of the static potential. These potentials have been calculated 
]46| , and in the limit of small distances appear to smoothly connect to the gluelump 
spectrum |T0|. In view of this the splitting of the lowest 1 and l +_ hybrid mesons, 


computed in the appropriate potentials, should be close to the gluon mass. However, 
these states have so far been elusive experimentally. 

The second possibility is a splitting between a hybrid and a non-hybrid state. Since 
spin effects are suppressed, let us consider QCD with scalar quarks as specified in 
Appendix |A|. We are interested in the correlator of the vector meson 


Vi(x) = Im (V (x)Ui<t>(x + 1 )) , 


J PC = r 


( 6 . 1 ) 


but now we wish to keep the bare scalar mass m 0 finite rather than integrating over the 
scalar fields analytically. The main difference to the static case concerns the subtraction 
of the contribution of the scalar fields. Away from the static limit they will not sit at 
the same point in an adjoint state, but rather form spatially extended singlet mesons, 
the lightest of which couples to 


S(x) = Re (0 t (x)0(a;)) , J PC = 0++. (6.2) 

Comparing the operator content of S, V, one may think oiV as a hybrid meson. The 
energy difference between vector and scalar mesons is accounted for by three different 
contributions: i) the total spin J = + J g , consisting of angular momentum of the 

scalars and the gluon spin; excitation of the it) scalar and in) gluon fields in higher 
quantum states. In the limit mo —> oo the scalars become static and hence their angular 
momentum is switched off. Furthermore, they are “quenched” to be external fields and 
cannot be excited into higher quantum states. The remaining contributions to the mass 
splitting are the gluon spin and excitations of the gluon field. In this case the vector 
meson becomes spin-exotic, i.e. its quantum numbers cannot be accounted for without 
gluon degrees of freedom. Thus, the mass difference between vector and scalar meson 
in the static limit is probing a gluonic excitation with the quantum number of the 
gluon. This suggests that the gluon mass should be approached arbitrarily accurately 
by the limit 

mA = lim \My — Ms]. (6.3) 

mo—too 
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If the limit exists, the last equation may even serve as an alternative definition of the 
gluon mass. However, formally it is not clear whether a smooth limit exists. Calculating 
the scalar correlator in the static limit in the continuum using Eqs. ( |A|1.3| ),( |A|1.5|) , one 
obtains 

(S(x)S(y)) g%(x,y) (Tr (U 0 (x,y)U%(x,y))) , (6.4) 


which has the same form on the lattice. The temporal Wilson lines combine to a 
unit matrix, hence the desired cancellation of their renormalization with that of the 
V correlator does not work and, strictly speaking, the formal static limit of My — Ms 
does not have a finite continuum limit. The reason is that static sources in a singlet 
state annihilate when they are at the same point. The only way for them to exist at 
zero separation is to combine their colour into an ( N 2 — l)-plet state, which we have 
used in the previous sections to calculate the static limit. However, this problem is 
absent for arbitrarily large but finite mo, and one may ask whether the value for the 
gluon mass can be smoothly approached by My — Ms for arbitrarily large m 0 . As we 
shall see in the next section, this is indeed the case in 2+1 dimensions. 


6.1. Heavy mesons for SU(2) in 2+1 dimensions 


We now turn to a test of this proposition in the SU(2) lattice gauge theory in 2+1 
dimensions. The scalar QCD action Eq. ( |A|1.1| ) can then be rewritten in terms of 
scalar 2x2 matrix helds $(x) and discretized as given in Appendix |B|. Results of 
detailed simulations for varying scalar mass mo are shown in Fig. |6|, which displays 
the scalar and vector mesons, together with the lightest scalar glueball and its first 
excitation. 

Note that scalar glueballs and mesons are indistinguishable by quantum numbers. 
To tell one from another a mixing analysis following |22], ^3] has to be performed. Quite 
generally, the glueballs and scalar mesons show very little mixing even when they are 
close in mass, and thus are easy to identify. An example is shown in Fig. |6| (right), 
where the matrix elements between the 0 ++ mass eigenstates + n and the operators used 
in the simulation Oi are displayed for mo = 1-0 g 2 . Here the operators are “smeared” 
versions of pure gauge and scalar operators, O t e {G = Tr F?-, S = 0V}- (For details 
of the mixing analysis in a Higgs model, cf. [221). Beginning with the lowest one, the 


three 0 ++ states are easily identified as glueball, scalar bound state and glueball, based 
on their operator content, and the situation is similar for the other parameter values. 

The dotted error bands in Fig. [6] (left) give the location of the lightest scalar glueballs 
in the pure gauge theory |3T|. It has already been reported for various scalar gauge 


models in 2+1 dimensions that the glueball spectrum deviates only at the percent level 
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from that in the pure gauge theory |22], [23], ph7[1 . Indeed the figure shows that already 
for bare masses m 0 > 2 g 2 the glueball spectrum assumes its pure gauge values to high 
numerical accuracy, indicating that the scalars have largely decoupled as dynamical 
fields. 


V30.5 

0 
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V20.5 

0 
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V 10 .5 

"O.o 0.5 1.0 1.5 2.0 2.5 3.0 

/ 2 0 
m 0 /g 

Figure 6: Left: The lowest states in 3d scalar QCD, simulated at (3 — 18 on 48 3 lattices. 
M g denotes scalar glueballs, M S y scalar and vector mesons. 

Right: Overlap ('F i |C>j|0) between the three lowest states 4+ in the 0 ++ channel and 
blocked operators, Oj G {S, G}. 



m =1.0g 


Also shown is the mass splitting between vector and scalar mesons, which is only 
slightly diminishing over the range of mo studied here. This is in accordance with the 
fact that we are close to the pure gauge limit. Note however, that in three dimensions 
the Coulomb part of the static potential is logarithmic. In a logarithmic potential, the 
level splittings of bound states of scalars calculated from a non-relativistic Schrodinger 
equation are also independent of the constituent mass |48|], so that constancy of the 
splitting alone in this case cannot indicate the vicinity of the pure gauge limit. 

According to the prescription Eq. ( |6.3|) , extracting the gluon mass requires an ex¬ 
trapolation of the splitting to infinite scalar masses. However, for heavier scalar fields 
the numerical errors grow rapidly and an accurate determination of the mass difference 
is increasingly difficult. This is well known from heavy quark physics and can in prin¬ 
ciple be cured by NRQCD methods [[49]]. Taking the largest value of m 0 in Fig. [|one 


reads off tua = 0.37(6). This is in excellent agreement with the static limit calculated 
in Sec. |5.1|. We then conclude that in scalar QCD in 2+1 dimensions the gluon mass 
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is observable to good accuracy in the mass splitting between vector and scalar mesons, 
which smoothly connects to the static limit. 


6.2. Heavy quarkonia in QCD 


Similar to the pure gauge physics in Table |2|, a close resemblance is known between the 
SU(2) Higgs models in 2+1 and 3+1 dimensions. Both phase diagrams are analytically 
connected, with mass spectra of Higgs and W-bosons or bound states of scalars and 
glueballs in the Higgs and confinement regions, respectively (for a review and refer¬ 
ences, see [Q). Moreover, in the symmetric region the confinement and screening of 


fundamental charges as well as the mixing properties of the gluonic flux with meson 
states are the same in both models fl5l|], and one might expect the results of the pre¬ 


vious section to carry over. However, there is also a significant difference concerning 
the short range Coulomb part of the static potential, which is now ~ r -1 . While a 
calculation analogous to the one in Sec. |6T] is beyond the scope of the current paper, 
it is intriguing to speculate about the situation in QCD, which provides us with heavy 
quarks to probe the gluon dynamics. 

We then consider the mass splitting of heavy quarkonia instead of scalar mesons. 
Away from the static limit the situation is complicated by the quark spins, the total 
meson spin now being J = L (/ + J g + S q . Spin independent meson masses analogous to 
the case of heavy scalar fields (and static sources) are obtained by averaging over the 
quark spin multipletts. In the standard spectroscopic notation (based on the quark 
model) the total meson spin is written as J = L + S q . Hence L = L 5 + J g , which is 
now playing the role of J = + J g in the scalar case. Thus, based on a comparison 

of quantum numbers, the splitting My — Ms in the theory with heavy scalars might 
correspond to the well known spin averaged mass splitting A 1P _! 5 , as put forward in 
a preliminary account of this work 


Using current experimental numbers 


one 


finds A!p_! 5 (cc) = 418.5 MeV and A ip-is(bb) = 416 MeV. 

However, it is not clear that these numbers are tied to the gluon mass. The fact 
that they are practically identical and quark mass independent cannot be taken as an 
indication of the proximity of the static limit, but is also accounted for by the effec¬ 
tively logarithmic non-relativistic bound state potential at the length scales relevant for 
these quarkonia^ |48], [54|. Increasing the quark mass pushes the states into the ~ r -1 
Coulomb region of the potential, in which mass splittings scale linearly with the quark 
mass. Hence the splittings of orbitally excited quarkonia do not have a static limit. On 
the other hand, this statement concerns only states within the non-relativistic quark 


2 I thank W. Buchmiiller for pointing this out to me. 
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model. In QCD there should be additional states, the hybrids, in which gluonic exci¬ 
tations account for the same quantum numbers, and one would expect the static limit 
of the splitting to exist. In general mixing analyses are necessary in order to unam¬ 
biguously identify the nature of an observed state. A four-dimensional analogue of the 
calculation in Sec. [nl] could clarify this question. 

7. Conclusions 

It has been demonstrated that an unambiguous, non-perturbative definition of a gluon 
mass is possible in terms of an eigenvalue of the Kogut-Susskind Hamiltonian of lat¬ 
tice gauge theory, which can be interpreted as the energy of a gluon coupled to static 
sources. This eigenvalue can be computed in several ways. It dictates the asymptotic 
exponential decay of non-local pure gauge observables which, in a particular gauge, 
reduce to the gauge held propagator. Using the eigenvectors of the covariant Laplacian 
for its construction, these observables are free of Gribov copies and strictly positive, 
with effective masses being upper bounds to the ground state. In practice these oper¬ 
ators are numerically feasible in Higgs regimes, but not in confining dynamics, where 
their non-local nature results in an almost exclusive projection onto torelonic states. 

However, the same eigenvalue also governs the asymptotic exponential decay of ratios 
of static gluelumps, in which gluon fields are bound to adjoint sources. This establishes 
a non-perturbative relation between the gluon mass and manifestly gauge invariant 
observables. These can be easily computed using local operators, with well controlled 
infinite volume and continuum limits. Furthermore, in three dimensions the mass 
splitting of vector and scalar mesons approaches the same quantity in the static limit. 
In four dimensions this might be modified due to the different nature of the Coulomb 
potential. 

Viewing the lattice field theory as a statistical system, it is clear that this gluonic 
mass scale plays a fundamental role in colour dynamics: 1 jrriA is the largest correlation 
length in the pure gauge system, providing an infrared cut-off for virtual states as well 
as setting a scale for the screening of colour interactions. For thermal physics, in 
particular, this should allow a non-abelian definition of Debye screening in complete 
analogy to QED. Moreover, the fact that mn in three dimensions is finite implies a 
non-zero magnetic mass which we have computed for SU(2), and thus the screening of 
colour magnetic fields in a plasma. A physical interpretation of the zero temperature 
gluon mass requires more work to investigate the static limit of meson mass splittings. 

We then conclude that a general gauge invariant description of colour dynamics in 
terms of partonic degrees of freedom should be possible, with numerous interesting 
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questions to be answered. First, it would be very desirable to solve the projection 
problem of the non-local gluon propagators and numerically confirm the relations be¬ 
tween their exponential decay and the static mesons established here. An obvious 
generalization would be to apply the same techniques to quark propagators, which 
might offer an alternative approach to non-perturbative quark mass renormalization. 
Similar methods applied to Polyakov loops should also allow to give a gauge invariant 
meaning to the finite temperature static potential with sources in an octet state. Fi¬ 
nally, one may hope to get a new handle on some proposed confinement mechanisms. 
For example, a gluon getting massive dynamically is in acccord with what is expected 
in a picture of the vacuum as a dual superconductor. 


A. Continuum pure gauge theory with static sources 

Let us introduce a complex scalar Applet 4> a (x), ot = 1, ...AT, coupled to the gauge 
fields, with the action 

S^[A,4>] = j d 3 x |(D #1 0)*(x)D a ,0(x) +mo0*(x)0(x)} . (1.1) 

Together with the pure gauge action, the theory corresponds to QCD with scalar 
“quarks”. In the limit mo —> oo the terms \Di<p \ 2 are suppressed and decouple from 
the above action. In this case the scalar fields propagate in time only, describing 
static sources coupling to colour electric flux of the gauge fields. The scalar propagator 
satisfying (—DqD 0 + G(x, y ) = S(x — y) is known exactly: 

G(x, y) = g Q {x, j/)<J(x, y )U 0 (x, y), (1.2) 

where g 0 (x,y) is the free scalar propagator without background field and U 0 (x,y ) the 
temporal Wilson line from x to y. 

We are interested in observables of the type 0(x i, X 2 ) = <j>'(xi)M(xi,x 2 )(f)(x 2 ), where 
M(x 1 ,^ 2 ) is a pure gauge matrix operator containing Wilson lines and/or covariant 
derivatives. Doing the Gaussian integral over the scalars one obtains for the correlation 
function (ay = (xj,0),yj = (x, : ,f)) 

(0(x 1 ,x 2 )0^y 1 ,y 2 )) = Z~ l f DA 0[A] e" 5 ™^ (1.3) 

with 

Z = J DA Z 4 , e ~ s ™W, Z 4 , = [det(-D 2 + m 2 )]" 1 , 

0[A] = Z ^ 1 j Dm* 0{x l ,x 2 )0\y ll y 2 )e- s ^. (1.4) 
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In the limit m o —s► oo the scalar determinant becomes a constant and 0[A\ represents 
a pure gauge quantity, 

0[A] = g 0 (xi, yi)go(x 2 , 2 / 2 ) Tr [M(x 1 ,x 2 )U 0 (x 2 , y 2 )M\y 1 , y 2 )U 0 (y u aq)]. (1.5) 

A prominent example is the string of colour flux ending on scalar charges separated by 
r = |xi — x 2 |, which is described by the operator M(x\,x 2 ) = U(xi,x 2 ). I 11 the limit 
mo —s► 00 the correlator goes over into the Wilson loop, 

{0{x^x 2 )0{ yi ,y 2 )) g*(t)(W(r,t)), (1.6) 


decaying exponentially with the energy of the string in the presence of sources. 

In the same way as the energy of a flux tube, the gluon energy can be probed by 
coupling it to scalar sources. In making contact to the lattice strong coupling gluon 
wave function Eq. ( |2.10|) , we consider the simplest gauge invariant operator with J PC 
quantum numbers of the gluon, 

Vi(x) = Im (<£+( x)Di<f>(x )) , J PC = 1~~. (1.7) 


The connection to the gluon propagator is established by integrating over the scalars 
analytically. Using Eqs. dnp.flnp, one finds for the correlator 


(Vi(x)V,(y)) -9o(*,y)(Tr \D i (x)U a (x,y)D i (y)Ul(x,y) 


( 1 . 8 ) 


Since the expression is manifestly gauge invariant, it may be evaluated in any gauge. 
Choosing the temporal gauge, A 0 (x) = 0, we have U 0 (x,y ) = 1 and D i (x)U 0 (x,y ) —> 
— 2igAi(x), The right hand side of Eq. (|1.8| ) is thus identical to the gluon propagator 
in temporal gauge times a free scalar propagator, 


| Ao=0 = Ag 2 gl(x,y) (Tr A i (x)A i ( 2 /))| y4o=0 . (1.9) 

Note that the pure gauge operator in Eq. (|1.8|) is simply the continuum version of 
Eq. (|3.6| ). It is then clear that its exponential decay is determined by the spectrum of 
the continuum Hamiltonian in temporal gauge. However, in the continuum, D t (x) is 
an operator defined at a point, which allows its correlator to be written down without 
the introduction of non-local functions of the gauge field. 


B. Lattice actions 

In this appendix we give the lattice action and parameters used for the simulations in 
the various parts of this paper. 
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SU(2) Higgs model, Sec. |4.1| : The lattice action for the fundamental representation 
Higgs model in 2+1 dimensions with 2x2 matrix field <f>(;c), cf. Eq. (ED, may be 
defined as 


S{u,i] = + 

1 rl i 2 ' 

H—Tr + Pr -Tr — 1 


( 2 . 1 ) 


The relation between the parameters of the continuum action Eq. (|4.1|) and those of 


the lattice action is at the two loop level floo 
4 


P = — 


Pr = x- -, 


y = 



_ 3 _ x | + (1 + Ax) 

8 \Ph P 32n [ + X) 


167T 2 


51 N 

-h 9x — 12x 2 

16 


I 11 + + C) + 5.0+ 5.2+ 


( 2 . 2 ) 


with the numerical constants £ = 3.17591 and ( = 0.09. Due to superrenormalizability, 
Eq. ( p.2|) are exact up to discretization errors. The continuum limit is at one point in 
the lattice phase diagram, Pq —> oo ,Ph = 1/3 , Pr = 0. 

Yang-Mills limit, Secs. |4.2| , |5.1| : The pure gauge action is seen to be a limit of 
the above obtained by taking Ph = Pr = 0. From Eq. Q2.2| ) it is evident that this 
corresponds to sending the continuum parameters x, y —> oo. 

Scalar QCD, Sec. |6.1| : Finally, the Gaussian scalar action simulated in Sec. 5+ is 
obtained by taking Pr = 0. 

The algorithm used to perform the Monte Carlo simulation using the action in 
Eq. (^T|) is the same as in |22], |23|. The gauge variables are updated by a combi¬ 
nation of heatbath and over-relaxation steps according to |56], [57|], while the scalar 


degrees of freedom are updated combining heatbath and reflection steps as described 
in 0. The ratio of the different updating steps is suitably tuned such as to minimize 
autocorrelations. In the simulations we gathered typically between 5 000 and 20 000 
measurements taken after such combinations of updating sweeps. 

In order to increase the overlap of operators with the low energy states to be mea¬ 
sured, standard “smearing” or “blocking” techniques have been applied to gauge and 
scalar fields. The correlation matrix between operators of different smearing levels is 























measured and diagonalized by a variational calculation allowing to extract ground and 
excited states in a given quantum number channel. For details and references see ^2 . 
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